Optimal Asymmetric Quantum Cloning 
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While the no-cloning theorem, which forbids the perfect copying of quantum states, is well-known 
as one of the defining features of quantum mechanics, the question of how well the theory allows a 
state to be cloned is yet to be completely solved. In this paper, rigorous solutions to the problem 
of M — ► N asymmetric cloning of qudits are obtained in a number of interesting cases. The central 
result is the solution to the 1 — > N universal asymmetric qudit cloning problem for which the exact 
trade-off in the fidelities of the clones for every N and d is derived. Analogous results are proven 
for qubits when M = N — 1. We also consider state-dependent 1 — > N qubit cloning, providing a 
general parametrization in terms of a Heisenberg star Hamiltonian. In all instances, we determine 
the feasibility of implementing the cloning economically, i.e., without an ancilla, and determine the 
dimension of the ancilla when an economic implementation is not possible. 



I. INTRODUCTION 

Quantum information constitutes a feasibility study of 
protocols subject to the physical laws of quantum me- 
chanics. When they cannot be realized, it quantifies and 
determines the best possible result. Often these consid- 
erations are made under restrictions such as those of lo- 
cal operations and classical communication. However, 
even in their absence, there is the fundamental restriction 
that any protocol must be implemented by a completely 
positive (CP) map. The quintessential instance of this 
problem is quantum cloning; it is impossible to take one 
copy of an unknown quantum state and produce a second 
identical copy, unless that unknown state is known to be 
drawn from a set of orthogonal states [l|. 

The no-cloning theorem, as a cornerstone of quantum 
theory, is related to many fundamental ideas such as the 
principle of no-signaling j2j, uncertainty relations and 
state estimation [3j. Apart from the intrinsic theoretical 
interest, the no-cloning theorem has also found applica- 
tion in quantum cryptography Q where it enables detec- 
tion of attempts by an adversary to copy the information 
on a communication channel. While the no-cloning the- 
orem is now well established, the question of the extent 
to which an unknown quantum state can be copied has 
been the subject of intensive research, excellent reviews 
of which can be found in [1, @] . 

The optimal cloning of discrete quantum states began 
with the idea of the Buzck-Hillcry quantum cloning ma- 
chine which obtains two identical copies of a given un- 
known qubit state 0] • This has been extended to M —> N 
cloning [8( where, starting from M copies of the same un- 
known quantum state, the task is to produce N output 
copies of as high a quality as possible. While the original 
cloning machines were symmetric, in the sense that all 
output copies had the same fidelity, this has also been 
extended to asymmetric cloning [9]. In this latter task, 
not all copies need to have the same quality; some out- 
put clones can be made to have higher fidelities at the 
expense of others. This is an important feature since it is 



precisely these trade-offs which quantify a potential level 
of eavesdropping, or the relationship between different 
correlations in a given quantum state. 

The original cloning machines were designed to clone 
all input pure states of a given dimension d (referred to as 
"universal cloning"). With prior knowledge of the distri- 
bution of the states to be cloned, better strategies can be 
devised for this, state-dependent, cloning. There are sev- 
eral well-known instances, the most common of which we 
refer to as "equatorial cloning" for d = 2, the state 
of the qubit is known to be drawn from the set of states 
on the equator of the Bloch sphere, (|0) + e'* |1))/V2 
and any angle <j) € (0, 2tt) is equally likely. This gener- 
alizes to "phase-covariant" cloning: for d = 2, the state 
of the qubit is drawn from the set cos 8 |0) + sin6*e 1 ^ |1) 
with a probability distribution that is independent of the 
parameter <f>. 

Another important question is the economy of imple- 
mentation of a doner; if the optimal cloning machine 
can be implemented unitarily without any ancillary sys- 
tems, the doner is said to be economical, otherwise it 
is not [III, [l2| . The presence or absence of any ancilla 
systems alters the implementation of the doner exper- 
imentally, the economic doner being simpler to control 
and less sensitive to decoherence effects. 

In this paper, we apply the powerful formalism of the 
Choi- Jamiolkowski isomorphism to solve several interest- 
ing instances of the cloning problem. The central result 
constitutes a completion of the study started in (l3| : a 
determination of the optimal cloning fidelities for 1 — > N 
universal asymmetric cloning of qudits (arbitrary N and 
d). While special case results are known for A up to 4 
@, Ell , we proposed in [l3[ a general solution, but some 
caveats remained unresolved. Wc simplify the proof of 
[lH , and extend it to close the potential loopholes in the 
argument, making use of the Lieb-Mattis theorem [15] to 
prove that the maximum eigenvalue of the whole system 
is found in a particular subspace. 

The paper is organized as follows. In Section|TTl we ex- 
plain the formalism used to obtain bounds on the achiev- 
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able fidelities in state transformation tasks and state the 
necessary and sufficient condition for the maximum fi- 
delity to be realized by a completely positive map. We 
also explain the applicability of the formalism to the 
problem of cloning quantum states, stating the figure of 
merit used to define the success of the operation. In sec- 
tion [nil we study the 1 —¥ N qubit cloning problem, pro- 
viding a parametrization for a class of state-dependent 
doners, deriving solutions in the symmetric instance and 
studying the economy of its implementation. In the 
asymmetric case, it is proven that if N is odd, the opti- 
mal doner is economical. For even N, no more than one 
ancilla qubit is required. For equatorial qubit cloning, we 
solve the asymmetric case for N — 2, enabling a simple 
derivation of the monogamy of the CHSH Bell inequality. 

Section IIVI constitutes the main result of the paper, 
where the 1 — > N asymmetric universal qudit cloning 
problem is solved and its economic realization is stud- 
ied. Many of the technicalities of the proof, which el- 
evate it over the conjectures in [l3j], are expounded in 
the Appendix. For N > 2, we prove that the optimal 
machine is economical. For N = 2, the optimal machine 
is not economical, requiring one additional qudit 1 . Sub- 
sequently, in Sec. |Vj the applicability of the formalism 
to the M — > N cloning is explained. The penultimate 
section, Sec. IVI1 contains the trade-off relations in the 
optimal fidelities in the two instances where such a re- 
lation can be derived, namely 1 — > N qudit cloning and 
(N — 1) — >• AT qubit cloning. The paper concludes with a 
summary of the work and possible future extensions. 



II. THE CHOI-JAMIOLKOWSKI 
ISOMORPHISM 

We begin with an explanation of the main tool used 
in the solution of optimal cloning tasks, the well-known 
Choi-Jamiolkowski isomorphism. This formalism is used 
in general to find how well a particular state transforma- 
tion task can be achieve d by a quantum process, i.e., a 
completely positive map (l 61 ] . 

The scenario is that, given one of a set of L states 
(i = 1,...,L), we are required to perform a particular 
state transformation on it, without being told which of 
the states we have been given. The required transforma- 
tion may not be achievable exactly within the quantum 
formalism, but is best approximated within the theory 
by a completely positive, trace preserving map £ that 
transforms input state \ipi) into £ The success of 
the state transformation task is then measured by a fi- 
delity given by 

F=i^Tr(M£(hM))- 

i 



1 other than extreme cases which reduce to N = 1 



Here Mi are positive operators (Mi > 0) satisfying 
II -Mi II < 1 so that F is indeed a fidelity taking values 
between and I. If the fidelity takes value 1, we in- 
fer that the map has perfectly implemented the required 
state transformation for all the specified input states. As 
a simple example, consider being required to transform 
the states \tpi) into states \(pi), in which case we simply 
define Mi = \4>i) (<pi\. 

The fidelity can now be rewritten in a manner that 
yields definite upper bounds. This is accomplished by 
the isomorphism as follows. Since £ is a completely pos- 
itive map, its operation on a subsystem O of a bipartite 
entangled state (entangled between input system / and 
output O) is well defined. Let us denote the maximally 
entangled state of interest as 

d-l 
* n— 

where the basis states \(j>i) (i = 0, . . . , d— 1) span the sub- 
space occupied by the set of input states \ipi) (which have 
dimension d). The action of the composite map, com- 
posed of the identity on the input space and the desired 
map £ on the output space of the maximally entangled 
state gives us the output xio, 

li®£{\B))o =Xio- 

The condition that £ be trace preserving then implies 
that Tr(x/o) = (B\B) = 1. If the map is not trace pre- 
serving, this would instead yield Tr(x/o) < f , which will 
not alter our conclusions since we are interested in finding 
an upper bound to the fidelity of state transformation. 
The fidelity can now be evaluated as 

F = | Tr (xio(m ® Mi)) , 

i 

where the superscript T denotes transposition. If we de- 
fine a matrix R as 

fl=£l>iMl<®Mi, (I) 

i 

then we obtain 

F = Tr(Rxio) < Tr( X /o) 'max _ ^max (2) 

where r max is the maximum eigenvalue of matrix R. We 
thus arrive at the important conclusion that the max- 
imum fidelity achievable within quantum theory for a 
given state transformation task is bounded by the maxi- 
mum eigenvalue of a suitable matrix R. 

A few points are worthy of emphasis. The dimension d 
appearing in the equations above is the dimension of the 
subspacc of Hilbert space spanned by the set of states 
Even if the states 1^) are defined on M qubits, it 
is not necessarily the case that d = 2 M , since the states 
might exist on a smaller subspace. Such is indeed the 
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case when we are performing the M — > N cloning task 
as we shall see in a section |Vj For a continuous set of 
states, the sum appearing in the definition of R (Eq. JT])) 
transforms to an integral. The factor of 1 / L appearing in 
Eq. (fTJ) stems from the assumption that each of the input 
states \ipi) is equally likely. If this is not the case, these 
parameters can be adjusted based on a given probability 
distribution of the input states. 

The fact that the state transformation problem has 
been transformed into the problem of finding the max- 
imum eigenvalue of a matrix is interesting and highly 
useful. Even if the matrix R proves to be difficult to diag- 
onalize exactly, one can use many techniques (borrowed 
for instance from condensed matter physics) to bound 
the achievable fidelities. One example of such a tech- 
nique that has been used in finding bounds on ground 
state energies of condensed matter systems is to upper 
bound the norm of the matrix by the sum of the norms 
of the constituent terms (t7| . 

A. Achieving the Maximum Fidelity 

From the equations above, we see that the maximum 
eigenvector, \*ff), of the matrix R defines the optimal 
strategy if it can be realized. If |\&) is unique, it can 
be expressed as a pure bipartite state between the sub- 
systems I and O with a Schmidt decomposition 

d-l 

|M/)=X>n|<^)|A n ). 
n=Q 

When the Schmidt coefficients are given by j3% = 4, 
the state is maximally entangled across the input/output 
partition and can be implemented by a unitary U, defined 
as 

U\4> n ) = |A„) Vn. 

Here the relevant Hilbert spaces are extended as neces- 
sary so that they have the same size. In this instance, 
the optimal strategy is called economical, meaning that 
one does not require an ancilla for the operation to be 
implemented. In fact, even if the maximum eigenvector 
is not unique, provided there exists a superposition of 
the maximum eigenvectors that is maximally entangled, 
the optimal map can be implemented as a unitary and is 
therefore economical. 

More generally, if there exists a mixture of the max- 
imum eigenvectors of R, pr, such that TtoPr is maxi- 
mally mixed, then this can be implemented as a CP map 
or, equivalently, a unitary operator over a larger Hilbert 
space, in which case the operation is no longer economi- 
cal. The condition for implementation of the state trans- 
formation task optimally by a CP map is therefore 

That this condition is sufficient is seen easily by recogniz- 
ing that one can add an auxiliary Hilbert space to purify 



Pr. The overall pure state then defines a unitary as in the 
previous case although this operation is not economical. 
That this condition is also necessary is seen by writing 
the optimal CP map using the Kraus decomposition as 

£(p) = J2^ P Al 

i 

where J2i = l - We trien nnd tnat 

p R = lj ® £ (\B) (B\) = £ A° \B) (B\ IO A?. 

i 

And upon performing the partial trace over O, we obtain 
Tr PR = Tr f^^|fl> (B\ IO J , 

which can be reduced using the completeness relation for 
the Kraus operators J2i A\Ai = 1 to TroPR = ^1/. 

As a final remark we mention that the optimal strategy 
can also be implemented by teleporting the input state 
onto spin / of a resource state which could be either 
|\&) or the purification of pr. The different measurement 
results of teleportation can be corrected for by action 
on the output space if the state satisfies the required 
trace condition. In the case of cloning, this gives rise to 
telecloning protocols. 

B. Application to cloning quantum states 

The potentially powerful formalism described above is 
now used in the problem of the optimal cloning of quan- 
tum states. In the cloning process, we start with an 
unknown quantum state of Hilbert space dimension 
d, and aim to produce N copies of the state. It is known 
that this state is drawn from a possible set of states E 
with distribution f(ip) that is normalized as 

f /(#ty = l. 

Dividing the output space O into N qudits labeled I to 
N, our aim is to optimize the quality of the TV different 
copies to be produced. There are different figures of merit 
that can be applied to cloning (i.e. different definitions of 
the Mi). The simplest figure of merit is the global fidelity, 
for which A4i = \ipi) (ipi\® N . The solution to the cloning 
problem in terms of the global fidelity is known [l8j and 
will not concern us in the rest of the paper. Instead, 
we consider the single copy fidelity, in which we take the 
global output state pi...N, and assess the fidelity of a 
single copy on a given site n, F n = Tr(|V>i) (i>i\ n pi...N)- 
To this end, we define 

N 

M l ^^a n \ipi) (ipi\ n 

71=1 
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where imposing a n > and 5Z n=1 a n = 1 ensures that 
the A^j satisfy the required property ||A^j|| < 1 in ad- 
dition to Mi > 0. In particular, F = 1 can still only 
be achieved if the output state is \ipi)® N for all inputs i. 
For generality, we assign different weights a n to the dif- 
ferent copies, emphasizing a possible desire for different 
qualities of output, although a common desire is equal 
qualities, a n = 1/N. The latter case is called optimal 
symmetric cloning while the general scenario is the opti- 
mal asymmetric cloning problem from 1 to N copies. 

For the single copy fidelity, on which we henceforth 
concentrate exclusively, the matrix R is seen to be 

N 

R= / /W#|V>><<® (3) 

We are considering here cloning from a single copy of \4>) 
so that the problem is 1 — > N cloning. This means that 
the dimension of the input Hilbert space I is no larger 
than d, and the basis is simply \i) for i = . . . d — 1. 

III. STATE DEPENDENT CLONING OF 
QUBITS 

In this section, we illustrate the above formalism by 
restricting to the common case of qubits (d = 2) , param- 
eterizing the input state by 

|V)=cos^|0)+sin^|l) 

for some distribution function f(6,(f>). We impose two 
restrictions on the input distribution function f(6,<j>), 
namely: (i) f{6, <j>) is phase covariant, meaning it is in- 
dependent of (f> } i.e., f{0,4>) = f{9) and (ii) f(0,<j>) is 
symmetric about the equator of the Bloch sphere, i.e., 
f(0) = f{n-0). 

Using these assumptions 2 , the matrix R for the 1 — > N 
cloning of qubits can be written as 

R = i i + J2 anT ( xx - YY + ^^- zz ) ' 

n=l V ZL /0,n 

where X, Y and Z are the familiar Pauli matrices, and 

r = i/ •^ 6, ' < ^ sin2 dddd<t) - 

The parameter, T, varying between and \ provides 
an intuitive interpretation regarding the distribution of 
states over the Bloch sphere - the larger the value, the 



more tightly packed the states are around the equator. 
We would now like to find the maximum eigenvector and 
maximum eigenvalue of R. The problem can be recast 
by applying a rotation Y to R, and instead demanding 
the minimum eigenvector (ground state) of a new matrix 
R given as 

R = -^ + Y, a " r \ xx + YY + ^^ zz ) ■ 

„=1 V ZL /0,n 

This matrix R is familiar as the Hamiltonian for an 
anisotropic Heisenberg model on a star configuration and 
we can proceed to find its ground state. 

Before we do that, let us analyze the possibility of im- 
plementing the above optimal qubit cloning transforma- 
tion economically. To do so, we observe that R commutes 

with the operator X)n=o^™' so tnat tne svs t em divides 
into a series of excitation subspaces. We also note that 
[R, X®( Ar+1 '] = 0, which maps the k excitation subspace 
to the N+ 1 — k excitation subspace. These two observa- 
tions enable us to make some general statements about 
the entanglement of the ground state, and hence about 
the economy of cloning. To do this, consider the mini- 
mum eigenvector in an M excitation subspace, which can 
be split between the input and output systems, 

|*M> =a\0)\4> M )+P\l)\4>M-i)- 

Since the corresponding state in an N + 1 — M excitation 
subspace, 

|*jv+i-m> = a |1> X® N I^m) + P |0) X® N \<t> M -i) 
is also an eigenvector, one can construct the state 

^(|*M> + |*JV+1-M» 

which is also an eigenvector with the same eigenvalue. 
It is readily verified to be maximally entangled unless 
M = \N, \{N + I) or \N + 1. For odd N, consider 
the subspace with ^(N + 1) excitations. If l^fjv+i)^) is 
the unique minimum eigenvector in this subspace, then 
invariance under X®( N+1 *> proves that it must be maxi- 
mally entangled. On the other hand, if this state is not 
the unique minimum eigenvector, there always exists a 
state in the degenerate space which is maximally entan- 
gled. Hence, for odd N, the optimal doner is always 
economical, meaning that no ancillary system is required 
to implement the cloning transformation. For even TV, 
this is not necessarily true, but no more than one addi- 
tional qubit is required because one can increase N by 
I and then implement the optimal doner, imposing that 
one of the fidelities is \ (ajv+i = 0). 



2 The assumptions specified have a natural interpretation, 
but in practice one only needs the following four inte- 



grals to vanish: / f(0, 4>)d0d4> cos 6, J f(0, 4>)dHd<f>e'^ sm t 
f f(0, ^ddd^ sin(20) and / f(0, <t>)ded<t>e i2, t> sin 2 0. 



A. Symmetric Cloning 

The most commonly studied instance of cloning is 
where all the output copies are required to have the same 
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fidelity, so a T: 



R 



Thus, we have 




The permutation invariance on qubits 1 . . . N means that 

the operators such as (jfJ2n=l re duce into a simple 

direct sum structure 3 . In the present instance therefore, 
the eigenvalue that we require will always be taken from 
the fully symmetric subspace, meaning that the problem 
reduces to solving for the minimum eigenvalue of 



X,N+1 



N 



YiS 



Y,N+1 



Z,N+1\ 



where S x ' is the X-rotation matrix for a k dimensional 
space (see the Appendix for further information) and 
A = 1 ~^ r . The maximum fidelity that can be realized 



2r 



in the cloning transformation is thus given by 
i TA 



F 



— max i V 4 (* + l )( N ~ i) + A2 ( N ~ 2i ~ !) 2 - 

The choice of the best i depends on the value of T. For 
\ < r < |, the i should ideally be (N - l)/2. How- 
ever, since i must be an integer, it takes the value [N/2\ . 
Hence, in this regime of T, the fidelity is given by 



F 



N odd 



For < r < i, the term is maximized when i is either 
or N — 1. Both give the same fidelity, 



F = - 
1 2 



i -4r 

2N 



+ —\JAN + (N - 1) 2 A 2 . 



For classical inputs, when the subset of possible states 
is only |0) or |1), then it is clear that one should be able 
to achieve the maximum cloning fidelity 1 . This is indeed 
the case, because when T = 0, the maximum eigenstates 
of R are \0f (JV+1) and \lf (/ " ' 
a maximally entangled superposition 



and one can construct 
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71 



|oy 



ID 



proving that there exists a unitary that achieves the op- 
timal fidelity of cloning. 

If the state is selected uniformly over the surface 
of the Bloch sphere (r = |), R becomes the Hamiltonian 
of the isotropic Heisenberg model. The optimal fidelity 
of symmetric universal cloning, F = 2 ~^ 1 , is recovered. 



3 When asymmetry is present, but subsets of spins should have 
the same fidelity, those sets have that direct sum structure, and 
are replaced by a single spin operator of the correct dimension. 



B. Equatorial Cloning 

In the case of equatorial cloning, the set of states is 
restricted to the input distribution function f(6) with 
the specific choice of 9 = ir/2. In this situation, we clone 
states that are on the equator of the Bloch sphere, T = j. 
While we do not solve this problem for general asymme- 
try, the fully symmetric case 10] is now easily recovered 
by our previous considerations. 



F = 



3W+1 
4 



y/N(N+2) 
4N 



N odd 
TV even 



We can examine the possibility of achieving the optimal 
cloning fidelity for equatorial qubit cloning in an eco- 
nomical manner even in the general case of arbitrary 
asymmetries. Even without completely obtaining the 
solution, we can make the following statement about 
the entanglement of the telecloning state - any optimal 
pure state must be maximally entangled. This strong 
statement, showing that the transformation can be per- 
formed economically, follows from the observation that 
Z (R+ \l)Z = -(R + Hence, if \$ M ) is a mini- 
mum eigenvector in a particular excitation subspace, M, 
then Zq I'J'm) is a maximum eigenvector in the same ex- 
citation subspace. With the exception of eigenvectors of 
R with eigenvalue — \ (which give a fidelity smaller than 
a complete guess), these two states are orthogonal, mak- 
ing the computational basis the Schmidt basis, with equal 
Schmidt coefficients. Therefore, we can conclude that for 
arbitrary asymmetries and any N, the optimal equatorial 
qubit doner is economical. This generalizes the results 
in which show explicitly that for the symmetric case 
and N = 2, the equatorial doner is economical. 

Serving as a prelude to the derivation of monogamy re- 
lations from solutions to optimal cloning problems in Sec. 
IVI[ we make a brief aside to indicate how the monogamy 
relation of the CHSH Bell inequality [19( arises as a con- 
sequence of the consideration of 1 —> 2 asymmetric equa- 
torial cloning. The CHSH inequality is a two-party two- 
setting Bell inequality. Up to local unitaries, the operator 
measured in a CHSH test on a quantum system between 
parties and n can be expressed as 

So» = V2(XX - YY) 0n . 

This is due to the fact that (up to local unitaries), the 
individual settings can be chosen to lie in the X-Y plane 
by both parties. Moreover, for optimal violation for any 
entangled quantum state, the local settings should be 
chosen such that they are as far away from being com- 
muting as possible |20|. In this scenario, we would like to 
derive the optimal trade-off between the CHSH inequal- 
ity violation between Alice and Bob (each holding a qubit 
state) and that between Alice and Charlie, i.e., between 
Bqi and Bq2- Maximizing the cloning fidelity corresponds 
to finding the maximum eigenvalue, A, of 



4^2 



(aiBoi + a 2 B, 



02 ) 



6 



which implies that the average values of the two CHSH 
parameters for any state, (Bqi) an d (B02}, obey 

ai (B i) + a 2 (B 2) <4^2(A-i) 

where the maximum eigenvalue is given by A = 
a i + a 2^~ !)• The two asymmetry parameters must 
satisfy ax + <x 2 = 1, so we can choose them to satisfy 



ai 



(B i)k a 2 = (B 02 )k, 



for suitable k. This choice of the asymmetry parameters 
implies the well-known monogamy relation for the Bell- 
CHSH inequality within quantum theory (l9j 



(B 



(Bo 



< 



In general, when interpreting the state \^>) as a state 
used in telecloning, its ability to clone well (and hence the 
general performance of the corresponding unitary doner) 
is limited by the monogamy of its correlations. This could 
be derived for any of the state dependent doners. Doing 
so for the universal doner produced singlet monogamy 
QJ] (see Sec. ED). 



IV. 1 -> N QUDIT CLONING 

We shall now move to the main result of the paper, 
namely the solution to the general case of 1 — > N cloning 
where the input state lives in an arbitrary Hilbert space 
dimension d, and the output copies can have differing 
fidelities. For this qudit cloning problem, we concentrate 
on the specific case when f(tp) is uniformly distributed 
over the set of states according to the Haar measure i.e. 
we are equally likely to be trying to clone any one qudit 
pure state. To derive the matrix R for this, the 1 — » N 
universal asymmetric qudit cloning, problem, we write 
that = U |0) and integrate over U to obtain [l3[ 

N 

R = d dUU* |0) (0| U T ® a ^ 1°) (°\n u - 

J n=l 

The integration results in twirling [2l[ to yield 



R = 



N 



1 



1 ^ 

71=1 



\B) (B 



0,71 • 



where \B) 



is the Bell state as de- 



1 

Vd ■ 

fined previously. 

The main route towards solving this is to notice that 
any symmetric state |$) of N — 1 qudits (meaning the 
same state is returned when any pair of qudits is ex- 
changed) allows one to write down states {l^m) = 
l^)o m l^^ 1 )! 2 Njtm) s P an a subspace of R. The 
permutation invariance ensures that 



\B) (B 



0,m 



l + (d-l)*„ 



Within the subspace, the eigenvectors of the matrix R 
are described by 



(4) 



subject to the normalization condition 

/ N \ 2 N 

+(d-l)Y J Pl = d. (5) 



Here the coefficients {/3 n } are chosen such that they form 
a vector which is the maximum eigenvector of the matrix 



N N 

2J a„ \n) (m\ + {d - 1) yj On \n) (n\ 



n=l 



(6) 



By definition, then, this is the maximum eigenvector of 
the particular subspace, eigenvalue A. If it were the max- 
imum eigenvector of the entire matrix R, the cloning fi- 
delity would be given by F < ^tj- To prove that this 
eigenvector is indeed the maximum eigenvector, we in- 
voke a variant of the Lieb-Mattis Theorem Y\M . Indeed, 
if d = 2, no modification is required, and |15| instantly 
proves that we have the maximum eigenvector because, 
by the previous results, we are searching for the ground 
state of a Heisenberg star. The full proof for arbitrary d 
is deferred to the Appendix. 

It should be noted that the matrix in Eq. ([5]), which 
is derived from the matrix R and used to choose the 
particular set of {/3 n }, is independent of the choice of 
|$). Also note that since the matrix in Eq. (|6]) is real, 
the coefficients /?„ can be taken to be real, justifying the 
absence of conjugates in the normalisation condition. 



A. Possibility of economic universal qudit cloning 

Having proved, as in the Appendix, that the state Eq. 
(|4]) gives the maximum eigenvector of R and that the 
optimal fidelity is given by F < -jty, we need to check 
that this bound can be saturated by verifying that the 
maximum eigenvector corresponds to an implementable 
CP map. To show this, we need to demonstrate that 
there exists a superposition or mixture of the set of op- 
timal states such that the partial trace over the outputs 
returns a maximally mixed state on the input fSec. |H|. 

Eq. ((4]) offers plenty of freedom in selecting the sym- 
metric state |$). Hence, we will pick a particular sym- 
metric state 



1 d-l 



V ®(JV-1) 



Provided N > 2, tracing out qudits 1 to N of the corre- 
sponding |$) just leaves 

Po = J 1 . 
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i.e. |^) is maximally entangled, so the optimal map can 
be implemented, and that implementation is economical. 

On the other hand, for N = 2, one can examine an 
arbitrary |$) and prove that the maximally entangled 
state only occurs if /?i/?2 = i.e. one of the clones will 
be perfect, and the other has no information about the 
input qudit at all; this is the trivial N = 1 case. Other- 
wise, economical optimal cloning if impossible for N = 2. 
Nevertheless, one can readily verify that the mixed state 
formed by the maximum eigenvectors 

d-i 



has the required partial trace po — i 1, meaning that the 
optimal fidelity can be realized with only one ancilla of 
dimension d. The results of , where it was proven that 
universal economical cloning is impossible for the 1 — > 2 
case, are clearly a special case of this result. 



V. M N QUBIT CLONING 

The study of the optimal cloning machines can be ex- 
tended to the case when, instead of a single input copy 
of the state to be cloned, multiple copies are provided. 
The matrix R can be derived for this scenario in a similar 
manner to the previous sections. The difference here is 
that the dimension d! of the input space is not simply 
d M . Instead, d' is the dimension of the fully symmetric 
subspace of the M qudits. The matrix R is therefore 
defined as 



N 



N 



R = d' /M# (g) H) (< m £ On |V) W n ■ 



As before, the maximum eigenvalue of R upper bounds 
the achievable fidelity of the cloning transformation and, 
if there exists some superposition of the maximum eigen- 
vectors which is maximally bipartite entangled between 
the input qudits and the output ones, the maximum fi- 
delity can be achieved by a unitary operation. 

Here, we concentrate on the special case of universal 
M — > N qubit cloning. After performing the twirling 
operation and projecting the input spins onto the sym- 
metric subspace, we obtain 



R — £^ a n R n 



observe that the set of states 
1 M 

1 V \ n ) |$ M \ \<f> N - M ) 



\^x 



71=0 



x k =0} 



forms a closed subspace. Here |n)j represents the basis el- 
ements for the input symmetric subspace and x £ {0, 1}^ 
is an iV-bit string of weight M. The superposition above 
therefore denotes that a symmetric state is the 

state of the M output spins for which the bit value in 
the bit string is 1 (of which n are in the state |1) and 
M — n are |0)). ~ M ) is the symmetric state of the 
remaining N — M output spins with i qubits in the state 
To verify that these states form a closed subspace of 
R, note that 



Rn \ipx,: 



M+2 



x n 1 



Here the sum over y is over all bit strings yglOjlj^of 
weight M with the condition that y — x in all except two 
places, one of which is bit n, y n = 1. The relation 



1 



M + l-x-z 



then implies the normalization condition. 

Since the {IV^i}} span the space, we can write 



(7) 



The matrix R within this restricted subspace, R, is writ- 
ten, for bit strings x, y £ {0, 1}^ of weight M, as 



Rx: 



M+2 
17+2 



, M+l 
r M+2 

En:i„=0 a n<5y„ 



X = y 

y = M — 1 



So, at least in principle, the maximum eigenvalue of this 
matrix can be found. The Lieb-Mattis theorem is di- 
rectly applicable, so we prove that the subspace contains 
the maximum eigenvector for the entire R, revealing the 
upper bound on the optimal cloning fidelity. It still re- 
mains to be checked that the optimal fidelity can be im- 
plemented by a CP map, which is more readily achieved 
on a case-by-case basis. 



VI. FIDELITY TRADE-OFF RELATIONS 



with 



2R n = 1 



1 



M + 2 



X,M+1 



X,, 



■Sf< M+1 Y n 



C Z,M+1 7 



The fidelities of the output copies are then given by F n = 
(*| R n |*) with the final fidelity being F = J2 n a n F n . 

The proof strategy exactly mirrors that of the 1 — > 
N universal cloning task (see the Appendix). We first 



The formalism outlined so far reduces the task of find- 
ing the optimal universal qubit doner to that of diago- 
nalizing an x (^) matrix, although we have not ex- 
plicitly performed this diagonalization in the asymmetric 
case. Finding the maximal eigenvalue specifies a trade- 
off between the parameters {a n } and the qualities of the 
clones given by F n = (^| R n Evidently, specifying 
the solution to the asymmetric cloning problem is best 
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done in terms of trade-off relations between the various 
output fidelities. In this section, we derive the fidelity 
trade-off relations in two cases of interest to us, namely 
I — > N universal qudit cloning and (N — 1 ) — > N univer- 
sal qubit cloning. 

Whenever the task of finding the maximum eigenvector 
of R can be reduced to a problem on an iV-dimensional 
space (where N is the number of different fidelities we 
might have in the asymmetric cloning), we have N pa- 
rameters /?i in the state 

AT 

l*> = £>hfc) 

n=l 

where {(V'i)} span the space. The fidelities F n are simply 
a quadratic function of the j3i, as is the normalization 
condition ( V E'| V I / ) = I. We can therefore expect to elim- 
inate the Pi in the normalization condition and replace 
them with the F n . Intuitively, one can understand that 
when the dimension of the space is equal to the number 
of asymmetry parameters a n , it is possible to explore the 
entire subspace corresponding to the optimal solution by 
varying the different parameters a n . If the dimension of 
the subspace is greater, this is certainly not possible, and 
we can only follow a constrained path within the space, 
as specified by the maximum eigenvector. At that point, 
extracting the optimal trade-off relation analytically be- 
comes challenging. Fortunately, in the two particular 
cases specified previously, it is possible to derive these 
fidelity relationships by the method outlined above. It is 
important to note that neither case ever actually requires 
the determination of the maximum eigenvalue of R. 



where U and V are arbitrary d-dimensional unitary ro- 
tations. It is directly related to the cloning fidelity in 
1 — > N universal qudit cloning [22| , 

_ p , n d + 1 
n ~ d+l ■ 

The previous expression can now be rewritten in terms 
of the singlet fraction as 

N N 
n—1 n—1 

Eliminating the variables (3 n from the normalization con- 
dition, the trade-off relation in terms of the singlet frac- 
tions for the I — > N universal asymmetric qudit doner is 
obtained: 

E^^^ + i^+brfs^) ■ (8) 

n—1 \n— 1 / 

This relation describes the 'singlet monogamy' [l3| . 
in comparison to the tangle monogamy [23| or Bell 
monogamy [l9j |. as it serves to bound the correlations 
that are achievable between different sets of spins in any 
quantum state. It can be seen that this single equa- 
tion encapsulates the known results on I — > N universal 
asymmetric cloning, subsuming previous results on I — > 4 
cloning and much more besides. 

B. Trade-offs for the (N — 1) —¥ N qubit doner 



A. Trade-offs for the I — ► N qudit doner 

For I — > N universal qudit cloning, we derived the 
normalization condition of Eq. ([5]). This means that a 
(possibly non-optimal state) must satisfy the inequality 



N 



N 



\n=l / n=l 

The cloning fidelities are given by 



F„ 



In order to eliminate the /?„ in favor of the cloning fideli- 
ties F n , observe that 



N 



N 



Vd{d+l)F n ~d =(d-l + N)J20n- 

n—1 n—1 

For simplicity, the trade-off relations can be expressed in 
terms of the singlet fraction, po, n - The singlet fraction of 
a pair of states between and n is defined as 

Pa . n = max (B\ U ® Vp , n U^ ® \B) 



Aside from M = 1 universal cloning, another instance 
for which the fidelity trade-off relation can be derived is 
(N — I) — > N universal asymmetric qubit cloning. Let us 
replace the (3 X that appear in the solution to this problem 
in Eq. (JT]) with 7„, where, since the bit string x is of 
length N and weight N — 1, there is only one zero bit, 
which we denote as being at position n. For this problem, 
the normalization condition reads 



iV 



JY 



2 > 



5> + E 



while each of the cloning fidelities can be found to be 

,/2 



F n = 1 



In 



The elimination of the 7„ therefore leads to the following 
fidelity trade-off relation for the case of (N — 1 ) — » N 
universal asymmetric qubit cloning 

N I N \ 2 

jv-i<£>„- ^yr 7 ^ , (9) 

n=l \n=l / 

describing the entire set of achievable fidelities. 
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Finally, to verify that these fidelities can be achieved 
by a CP map, observe that an equal mixture of the two 
states |^o) and as in Eq. (J7|), yields the identity 

on the N — 1 input spins when the N output spins are 
traced out, 

|T> 1 ...^(|*o><*o| + |*i)<*i|) = ^l. 

This implies that this cloning transformation uses one 
ancilla qubit (i.e. is not economical) but gives the requi- 
site trace satisfying the necessary and sufficient condition 
for its implementation by a CP map. 

VII. CONCLUSIONS 

The efficacy with which a given quantum state can 
be cloned is an important problem in quantum informa- 
tion, having implications in diverse fields such as quan- 
tum cryptography, state estimation, and entanglement 
monogamy. In this paper, we tackled the problem using 
the Choi-Jamiolkowski isomorphism, and provided solu- 
tions in a number of interesting scenarios. The investi- 
gation of 1 — > N qubit cloning yielded a parametrization 
for a large class of state dependent doners including the 
equatorial and universal doners. If N is odd, the optimal 
doner is economical and when N is even, no more than 
one ancilla qubit is required. As an interesting aside, we 
derived the monogamy relation for the CHSH Bell in- 
equalities in terms of the fidelity trade-offs in the 1 to 
2 asymmetric equatorial qubit cloning problem. For the 
general 1 to N asymmetric equatorial cloning problem, 
we proved that for any N and any asymmetry, the doner 
is economical. 



The main result of the paper is the solution to the 
1 — > N universal asymmetric qudit cloning problem, for 
which we derived the exact trade-off relations in the fi- 
delities of the output clones, Eq. ©. For N > 2, the 
optimal doner is economical while for N = 2, it is not, 
requiring one additional ancillary qudit for its implemen- 
tation. Generalizing to M — > N cloning, we derived the 
(tight) trade-offs in the fidelities for the special case of 
(N — 1) to N universal cloning of qubits, Eq. (J9]), and 
provided some insight into why it is these two cases that 
yield such a straightforward monogamy relationship. 

One aspect that we have not touched upon here is the 
efficiency of implementation of the optimal doners. How- 
ever, states such as those of the 1 — > N doner, Eq. (Q}, 
have a concise description in the computational basis, 
which is sufficient to guarantee that they can be easily 
prepared, i.e. with circuits that scale as a polynomial in 
N [24]]. While such an argument says nothing about the 
economy of the implementation, as ancillas may be used 
in intermediate steps, economy can also be achieved. 

The formalism presented here, and, primarily, the tech- 
niques for proving optimality, can potentially be applied 
in many other scenarios. A natural generalization in- 
volves the cloning of mixed quantum states, a problem 
known as broadcasting [25[. There has also been a vast 
amount of literature on the problem of cloning for con- 
tinuous variable systems (2|| to which the methods may 
be directly applied. Finally, it is potentially interesting 
to consider the optimal cloning of quantum observables 
and entanglement rather than entire quantum states |27| . 
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Appendix A: Derivation of Optimal Universal Qudit 
Cloning 

In this Appendix, we provide the optimal solution to 
the 1 — > N universal qudit cloning problem. In partic- 
ular, we find the maximum eigenvector of the matrix R 
which is defined for this problem as 

1 d N 

n—l 

A major ingredient in the derivation of our results is 
the understanding of symmetries in a system, particu- 
larly those of rotational and permutation invariance. We 
will recap the bare minimum here in order to familiarize 
the reader with our notation, but for further details, the 
excellent treatment in [28[ can be relied upon. We start 
by defining the spin operators 

d-2 

S x - d = J2V(n + l)(d-n- l)(|n) (n + 1| + |n + 1) (n\) 

n=0 
d-2 

S Y > d = ^2y/(n + l){d-n- 1)H \n) (n + 1| + i \n + 1) (n\) 

71 = 

d-l 

S z < d = ]T(d-l-2n)|n)H, 
n=0 

where the as yet unused subscript will be used later to de- 
note which spin the operators act on. If d = 2, these op- 
erators will be simply written as X, Y and Z respectively, 
recovering the usual Pauli operators. In the 1 — > N 
cloning problem, we have one input spin (I) of dimension 
d' , and N output spins (1 . . . N) of dimension d. These 
will require us to specify the total spin operators 



d'.d 



X 



gX,* + £ S X„ 



N 



J 



N 



n=l 
N 



jd'.d _ Q Z,d! , Q Z,a 

z i 2.^1 « 



n=l 

Jx + Jy + Jz 



Note that we have an unusual combination of negative 
signs in these definitions. This is due to the symmetry 
structure not being the usual invariance under the uni- 
taries 1J® N+1 . However, note that under the action of 

d-l 

71=0 

on the input spin, this is transformed into the usual to- 
tal spin operators. For this reason, we instantly recover 
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the commutation relation [Jz,J 2 ] = 0. Since R simul- 
taneously commutes with both J 2 and Jz, it has two 
quantum numbers S (where 4S(S + 1) is an eigenvalue 
of J 2 ) and Mz (taking values —S to S in integer steps), 
providing a lot of information about its structure. 

We now use the permutation invariant states of M qu- 



dits of fixed excitation number k = . . . M(d — 1) which 
are denoted as l*^)- In the case of d — 2, these are the 
uniform superposition of all states of M qubits contain- 
ing k excitations. For a set of M d-dimensional spins, we 
note their relations to the spin operators, 



(E^W) = (M(d-l)-2k)\*ff) 

\m=l / 



\m— 1 / 

Under the action of the individual terms in the matrix i?, these states transform as 



IRWRI IR\ I tftiV— 1\ _ J \ B )o,n \$i )l,2...N 7 tn 71 711 

\B) (B| ,„ \B^ m |*, ) 12 ^ m - | , ]B)o n K -i }i ^ n n ^ m 



r 



This therefore implies that the set of states 
{\B) m 2 N _£ m } spans a subspace of R for each 

i, corresponding to a value of M z = \{d— 1){N — 1) — i. 
Within this subspace, the eigenvectors are given by a 
superposition 

N 

m-E^knM- 1 )^...^, (Al) 

n=l 

and the coefficients must satisfy the normalization con- 
dition 

/ n \ 2 N 

K>» +(d-i)E^ = rf - 

\n=l / n=l 

Moreover, the coefficients {/3 n } are chosen to be the max- 
imum eigenvector of the matrix 

N N 

E a n \n)(m\ + {d-l)J2 a n\n)(n\, (A2) 

n,m— 1 n— 1 

this being independent of i. The state Eq. (|A1[) formed 
using the coefficients /3 n from the maximum eigenvector 
of the above matrix, is the maximum eigenvector of R 
within a given subspace i. The cloning fidelity using the 
maximum eigenvector is then given in terms of the max- 
imum eigenvalue A as F = ^tj- Before we proceed to 
find the absolute maximum eigenvector of R over all sub- 
spaces we examine the known case of symmetric cloning. 



1. Symmetric Cloning 

For the symmetric doner where each of the output 
clones in the 1 — > N universal qudit doner have the same 
quality so that a„ = we can readily proceed to calcu- 
late the (3 n and F. All the /?„ are equal, and hence, by 
the normalization condition, we have 

pn = \j 1 N 2 + N(d-r) 

and the final fidelity is 

1 2(N-l) 

F = h -. 

N N(d+l) 

This is known to be the optimal fidelity of the symmetric 
universal qudit doner [5], so even without finding any 
other eigenvalues, we infer that these eigenvectors (given 
by Eq. (jAll) with equal j3 n ) must be the maximal ones. 
These states, which we denote by |^ ym ), will be used 
later in the proof. 

Let us make an important observation about these 
states. The eigenvectors l^ 5 " 11 ) within a particular Mz 
subspace can be chosen to have a positive amplitude on 
every basis state except for those special cases in which 
the qudit state on the input spin does not appear in 
the rest of the string. For these special cases, such as the 
state |31210), which does not have a |3) on the any of 
the sites except spin 0, the amplitude of the eigenvector 
|4'^ ym ) is 0. These special cases can be ignored because 
they are eigenvectors of R with eigenvalue l/(d+ 1) (for 
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any arbitrary set of {a n }), and are hence never optimal. 
It is important to note that not every choice of _1 ) 
has this property of having a positive amplitude on ev- 
ery basis state within the particular subspace but one 
can always choose the eigenvectors l^^" 1 ) to have this 
property. For example, choices of |$|) could be |11) or 
(|02) + |20))/\/2 for d > 3, and these do not have positive 
amplitude on every basis state of the relevant Mz, but 
there are superpositions which do. 



two components, Rq and R±, where i?o contains all the 
diagonal elements, and R\ is the remainder. In the 
computational basis, |a), we have (a\Ro\a) = e a and 
(a | Ri \b) = K a b. Note that K a b > 0. Assume that in a 
particular excitation subspace, Mz, we know the maxi- 
mum eigenvector, 

a 

with eigenvalue Em z ■ Hence, 



2. Asymmetric Cloning 

It is now our task to prove that the maximum eigenvec- 
tor given in Eq. (jAll) within this subspace of given 
i is always the maximum eigenvector of R, for arbitrary 
asymmetries. To do so, we are going to make use of a 
modified Lieb-Mattis theorem [15[. The idea behind the 
proof is the following: we will prove that the maximum 
eigenvector within a given Mz subspace (making no re- 
striction on S) has non-negative coefficients on all the 
basis states. Since (apart from the discounted special 
cases) l^f m ) has positive amplitudes on all these basis 
states, this implies that 



E 

6 



Kb a fb = (Em z — e a )fa- 



Any other state must have a smaller expectation value of 
R, unless it is also a maximum eigenvector. Let us first 
try a state \a) as the ansatz. This reveals e a < Em z - 
Hence, we can take the modulus of the above equation, 



J2 K baft 



(E Mz - e a )\f a \. 



Next, we try the ansatz state 



> 0. 



^)=£i/.i 



Therefore, the 



\S&i) must be drawn from the subspace 

{\ B )o,m \®i ~ X )i,2...Njtm}> and hence must be the maxi " 
mum eigenvector that we solve for in Eq. (|A2j) . 

This would successfully establish that for the subspaces 
M z = -±(d-l)(A r -l),...,±(d-l)(AT-l),weknowthe 
maximum eigenvector. The maximum eigenvector in the 
other subspaces cannot have a larger eigenvalue due to 
the following argument. If it were true that these other 
subspaces had an eigenvector with a larger eigenvalue, 
then this eigenvector must be drawn from a subspace of 
S > |(d — 1)(JV — 1). However, any such eigenvalue 
is degenerate in Mz so it would also be present in all 
other Mz subspaces from — S to S, in particular in the 
ones in which we have already found a different maximum 
eigenvector. So, by contradiction, this establishes that we 
have the maximum eigenvalue, F = (1 + X)/(d + 1), and 
that it is (at least) (d — l)(N — 1) + 1 fold degenerate. 

This leaves us with the task of proving that, within 
a given Mz subspace, the coefficients of the maximum 
eigenvector are non-negative. To do so, divide R into 



This requires that there is at least one non-zero f a such 
that 



b 



K ba \f b \ < (E Mz -e a )\f a \ 



but since J2b K ba\fb\ > \ J2b K bafb\, this can only be sat- 
isfied with equality for every a, meaning that the coeffi- 
cients of the maximum eigenvector in each Mz subspace 
satisfy 

fa>0, 

up to a global phase factor. This completes the proof that 
the maximum eigenvector within a given Mz subspace 
has non-negative coefficients on all the basis states. By 
the arguments outlined previously, this establishes that 
the state Eq. (|A1[) is the maximum eigenvector of R, 
the cloning fidelity being given in terms of the maximum 
eigenvalue A by F = -g±y. 



